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ABSTRACT 

We present our early results of the mapping of the nucleus of the starburst 
galaxy NGC 253 and its immediate surroundings using the Infrared Spectro- 



(N ! graph ^ onboard the Spitzer Space Telescope. The map is centered on the nucleus 

I I of the galaxy and spans the inner 800 x 688 pc^. The spatial distribution of the 

^ I [Ne III] line at 15.55 yum and the polycyclic aromatic hydrocarbons feature at 

Tjj- I 11.3 /im peaks at the center while the purely rotationnal transition of molecular 

^ I hydrogen at 17.03 /xm is strong over several slit positions. We perform a brief 

investigation of the implications of these measurement on the properties of the 

_ star formation in this region using theories developped to explain the deficiency 

^ [ of massive stars in starbursts. 



Subject headings: Molecular hydrogen, Galaxies, Starburst, Ionization, Massive 
stars, Starburst Galaxies 



1. Introduction 

NGC 253 is a well studied nearby starburst galaxy located at a distance of 2.5 Mpc 
(Pence 1980) and is mainly powered by an episode of central star formation. This episode 
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is generally believed to be caused by material brought to the center of the galaxy by a bar 
which is clearly visible in the 2MASS image (Jarrett et al. 2003). The presence of a bar 
from both morphological and kinematical evidence, and its role as starburst trigger has been 
extensively discussed in the literature (Das, Anantharamaiah, & Yun 2001 and references 
therein, Scoville et al. 1985, Engelbracht et al. 1998). 

Four luminous super star clusters were discovered in the central region by Watson et 
al. (1996) with the Hubble Space Telescope. They derived a bolometric luminosity for the 
brightest cluster of 1.3 x IO^Lq which corresponds to 10% of the luminosity of within a region 
of a radius of 15" at the nucleus. Also, observations of Ha emission and earlier Einstein and 
ROSAT x-ray data (see Ptak et al. 1997, and references therein) revealed a starburst- driven 
wind, emanating from the nucleus along the minor axis of the galaxy. This wind was well 
delineated by Chandra (Strickland, Heckman, Weaver, & Dahlem 2000). Engelbracht et 
al. (1998) studied NGC 253 using near-infrared and mid-infrared lines and concluded that 
the properties of the Initial Mass Function (IMF) in the starbursting region is similar to a 
Miller-Scalo IMF that has a deficiency in low mass stars. Also, Boeker, Krabbe, & Storey 
(1998) and Keto et al. (1999) presented maps of [Ne ll] 12.81 ^m. and PAH emission. 

Observations in the mid-infrared from space by the Infrared Space Observatory (ISO) 
were used to identify fine structure lines, molecular hydrogen lines and emission from poly- 
cyclic aromatic hydrocarbons (PAH) (Sturm et al. 2000; Rigopoulou et al. 2002; Forster 
Schreiber et al. 2003; Verma et al. 2003). These authors mainly studied the nuclear region 
of the galaxy with the Short Wavelength Spectrograph (SWS) and ISOCAM. Rigopoulou et 
al. (2002) detected several mid-infrared H2 lines in the SWS spectrum of NGC 253. They 
derived an H2 (0-0) S(2)/S(l) temperature of 380 K with aperture sizes of 14"x20" and 
14"x27" for S(2) and S(l) respectively. Thornley et al. (2000) measured the [Ne lll]/[Ne ll] 
line ratio using SWS and found a value of 0.06. Since the Ne atom and Ne"*" ion have an ioni- 
sation potential of 21.56 and 40.95 eV respectively, this ratio traces the massive star content 
because it is very sensitive to the spectral shape of the UV radiation field. Thornley et al. 
(2000) addressed specifically this issue by modeling the neon ratio as a function of the upper 
mass cutoff. They used evolutionnary synthesis code an photoionisation code to predict the 
neon line ratio around a cluster of stars. Their results show that a [Ne lll]/[Ne ll] line ratio 
of 0.06 is consistent with an upper mass cutoff around 25 Mq. However, they found that for 
short-lived bursts (1 Myr), the measured neon line ratio is consistent with the formation of 
massive (50—100 Mq) stars. Verma et al. (2003) also analyzed the SWS fine-structure line 
mainly in the context of gas-phase abundances. Their neon line fluxes indicate a slightly 
higher [Ne IIl]/[Ne II] line ratio of 0.08. They found a correlation between the neon ratio 
and metallicity among the SWS starburst sample but excitation for a given metallicity was 
relatively low compared to local H II regions. This problem with the deficiency of massive 
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stars in starbursts was also found by Puxley et al. (1989) and Doyon, Joseph, & Wright 
(1994). 

In this letter, we present spectral maps of [Ne ill], H2 (0-0) S(l) and the PAH feature 
at 11.3 /im. The sensitivity and small aperture size of the IRS gives us the opportunity to 
do spatial mapping over small scales of the mid-infrared properties of the central region of 
NGC 253. A total of 45 positions were observed around the bright starburst seen at the 
center of the galaxy. These observations allowed us to build the map of the [Ne lll]/[Ne ll] 
line ratio where a spatial gradient is detected. The observations and data reduction are 
presented in § 2 and the maps are presented in § 3. The spatial differences in the different 
maps are briefly discussed in § 4. 



2. Observations and data reduction 

The observations were carried out during the science verification phase of the Infrared 
Spectrograph (IRS) onboard the Spitzer Space Telescope (Werner et al. 2004) using the 
Short-High (SH) module (Houck et al. 2004). The high sensivity and narrow slit sizes of 
the IRS were used to map of the central region of NGC 253. The spectra ranges from 10 
to 20 /im with a resolution R k. 600. Although the observations of NGC 253 were non- 
contiguous on the sky, we have produced a "sparsely sampled map" to trace the changes in 
warm gas and dust properties across the inner starburst. The map is centered at a position 
of RA = 00M7™33.2^ DEC = -25°17'19" ( J2000) and the step between each pointing position 
is 15" along the length of the slit and 9" along its width. 

The map, which consists of 45 different pointings, is centered on the nucleus of the 
galaxy and spans an area of 70'.'8 x 60'.'8 at a position angle of 271 degrees. This corresponds 
to the central 800 x 688 pc^ of the galaxy at the assumed distance of 2.5 Mpc. Each slit 
covers an area of 11'.'8 x 5'.'3 (133 x 60 pc). The exposure time for each slit position was 88 
seconds. 

The basic data reduction was performed by the standard IRS pipeline version 9.5 at the 
Spitzer Science Center (SSC). The pipeline removes detector artifacts and cosmic ray signa- 
tures and applies the dark and flatfield corrections (Spitzer Observer's Manual^, Chapter 7). 
Standard full-slit spectral extraction and flux calibration were performed as described in the 
SOM. The uncertainty in the line fluxes is dominated by the absolute flux calibration. At 
the present state, these uncertainties are on the order of 30%. Uncertainties on the absolute 
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fluxes have little impact on the uncertainties of the relative fluxes. Line ratios are mostly 
affected by the determination of the continuum near the emission line. In our dataset, this 
source of error is below 5%. 
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Fig. 1. — The upper line is the spectrum of the central position of the IRS spectral map 
which corresponds to the location of the super star clusters of Watson et al. (1996). One 
of the spectrum off-center (lower line) is shown for comparison. The units on this spectrum 
are scaled up by a factor 40. 
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The 12.8 /iin [Ne ll], 15.55 /xm [Ne ill] and the 17.03 fim H2 (0-0) S(l) lines are detected 
with a signal to noise (S/N) higher than 3, our detection limit, on all the spectra. Line 
fluxes were measured by fitting a gaussian profile after substracting a second order baseline 
fit. Table 1 presents the values of the line fluxes for the central position of the map, which 
corresponds to the nuclear region of the galaxy. Figure 1 shows the spectrum at the central 
position and a typical off-center spectrum. The values measured for the central spectrum 
are shown in Table 1. 
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Table 1. Line fluxes at the nuclear position of NGC 253 



Line Rest wavelength (^m) Total Flux(10~^° W 



PAH 


11.3 


PAH 


12.0 


H2 0-0 8(2) 


12.27 


PAH 


12.7 


[Ne ii] 


12.81 


[Ne III] 


15.55 


PAH 


16.5 


H2 0-0 8(1) 


17.03 


[S in] 


18.71 



228±68 
20±6 
5±1 

159±48 

147±44 
12±3 
19±6 
6±2 

40±12 



The spectral maps are shown on Figures 2 to 5. The total value of the flux in a given 
line is shown as the value for the whole slit. This way, each staring observation becomes a 
pixel and the result is a 5 x 9 sparse map. All the images are overlaid on the contours of the 
2MASS K-band image. 



3. Analysis 

The [Ne ill] map (Figure 2) shows that the peak emission from this line is located at the 
center of the galaxy. Further out, the absolute line flux drops quickly. The same behavior 
is also seen with the PAH 11.3 fira feature (Figure 3). However, the intensity of the H2 
S(l) line is spatially distributed quite differently. We can see from Figure 4 that the overall 
intensity of the H2 (0-0) S(l) line follows the shape of the major axis of the emission seen 
in the 2MASS image. 
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Fig. 2. — Spectral map of Ne ill overlaid on the K band contours from the 2MASS image. 
The orientation of this map and all the subsequent maps is equatorial. The position angle 
of the major axis seen in the 2MASS image is 51° (Pence 1981). 
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Fig. 3. — Spectral map of the PAH feature at 11.3 fun. 
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Fig. 4. — Spectral map of H2 (0-0) S(l) over the 2MASS K band contours. 
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Figure 5 shows that the [Ne lll]/[Ne ll] hne ratio increases by a factor 4-6 from the center 
of the map toward the edges. This increase is higher than our uncertainties on the relative 
fluxes (see section 2). The line ratio at the central pixel is 0.08, in agreement with what was 
found with ISO-SWS (Verma et al. 2003). The ratio peaks at 0.46 on the south-east edge of 
the map. 
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Fig. 5. — Spectral map of the [Ne lll]/[Ne ll]line ratio over the 2MASS K band contours. 



-14- 

4. Discussion 

4.1. Line fluxes 

The [Ne ill] map shows that the strong neon emission is concentrated in the nucleus and 
decreases sharply outside that region. This behavior is also seen in the [Ne ll] line (not shown 
here). It is likely that the super star clusters found by Watson et al. (1996) are responsible 
for the major fraction of the high [Ne III] flux. The PAH features also follow this trend. 
These are thought to be excited in star-forming environments where PAH emission arises 
predominantly in photodissociation regions at the interface between ionized and molecular 
gas (Verstraete et al. 1996; Cesarsky et al. 1996; Crete et al. 1999). 

The map of Figure 4 shows a clear correlation between the H2 (0-0) S(l) line and the 
major axis of the emission seen in the 2MASS image. This correlation is also seen in the 
H2 (0-0) S(2) line (not shown here) and indicates that an excitation source is located along 
the major axis and rises the H2 temperature enough to sustain an elevated emission of the 
17.03 /xm and the 12.27 /xm lines. The purely rotational transitions of molecular hydrogen can 
be triggered either by radiation or shocks (Burton, Hollenbach, & Tielens 1992). In dense 
(n > 10^ cm^"^) photodissociation regions, UV photons from massive stars can produce 
enough energy to raise the temperature of the H2 gas so that collisional excitation lines are 
seen. The molecular gas can also be excited by x-rays (Rigopoulou et al. 2002) or cosmic 
rays (Bradford et al. 2003). The fact that the H2 S(l) line and the PAH 11.3 /xm line don't 
trace each other may be indicative that different mechanism maybe at work at the central 
position and along the major axis. 



4.2. Spatial variation of the ionization field 

The observed increase of the [Ne lll]/[Ne ll] line ratio (Figure 5) implies a harder ion- 
ization fleld outside of the nuclear region of NGC 253. To flrst order, this means that the 
radiation fleld outside the nucleus is harder than away from the nucleus. In this section, we 
breifly review the various physical conditions that may lead to an increase of the value of 
the neon line ratio. We discuss the impact of a metallicity gradient, a modifled IMF, and an 
age gradient on the physical properties of the nuclear region of NGC 253. 

Martin-Hernandez et al. (2002) showed that decreasing metallicity leads to a hardening 
of the ionization fleld since it affects the line blanketing and the characteristics of the stellar 
wind. However, a large gradient is needed to account for the observed change. The data 
of Martin-Hernandez et al. (2002) show that in order to triple the neon line ratio, a change 
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in metallicity close to 1 Z0 is required. The maps presented here only cover a radius of 
0.5 kpc so the implied metallicity gradient needed to reproduce the data is on the order of 
2 Zokpc-i. 

Several authors (Puxley et al. 1989; Doyon, Joseph, & Wright 1994; Engelbracht et al. 
1998; Doherty, Puxley, Lumsden, & Doyon 1995; Achtermann & Lacy 1995; Beck, Kelly, & 
Lacy 1997) proposed a truncated or steeper IMF to explain the deficiency of massive stars 
in starburst galaxies. This would imply that the properties of the IMF would change from 
the nucleus to the outer disk. However, the gradient would imply variations in the IMF only 
if the stellar populations in the different regions were nearly the same ages. 

The analysis of Thornley et al. (2000) showed that ageing can cause a softening of the 
ionization field where a starburst occurs. They pointed out that the observed properties of 
short-lived starbursts are very sensitive to the massive star content within the first 10 million 
years of a star formation episode. This means that the hardness of the radiation field, and 
thus the [Ne lll]/[Ne ll] ratio, should decay rapidly if the starburst episode responsible for the 
ionisation field is short-lived. Using this model, our results would indicate that the outside 
regions are, on average, younger than the nuclear region, consistent with star formation 
propagating outward. 



5. Conclusions 

We present spectral line maps of [Ne ill], H2 (0-0) S(l) , PAH 11.3 fira feature as well 
as the [Ne lll]/[Ne ll] ratio and the warm molecular hydrogen temperature of the central 
region of NGC 253 and its immediate surroundings. The strength of the H2 lines follows the 
morphology of the major axis while the [Ne III] line and the PAH 11.3 /im feature peak at 
the location of the nucleus and show no strong correlation with the major axis, suggesting 
different physical mechanism in and outside the nuclear region for the PAHs and neon on 
one side and the molecular hydrogen on the other. 

A significant gradient in the [Ne lll]/[Ne ll] ratio is detected. The ratio rises toward the 
edges of the map with the lowest value being at the location of the nucleus. To first order, 
this indicates that the radiation field is harder outside the nucleus which in turn suggests 
that there is a relative deficiency of massive stars in the nucleus of the galaxy. With our 
dataset, we were able to briefiy explore the variation of the ionisation field within a single 
starburst galaxy and compare it with models that are usually compared to the ionisation 
fields of a sample of starburst galaxies. The spectral-mapping capability of the IRS allows 
detailed studies of the local environment of powerful starbursts and provide unprecedented 
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insight into the star formation process in starbursts. We presented here a small sample of our 
dataset. A more detailled study on the impact of our observations on the physical properties 
of the central region of NGC 253 and the stellar population is in preparation. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. Support for this work was provided by NASA through Contract Number 
1257184 issued by JPL/Cahech. 

Wc would alos like to thank the anonymous referee whose suggestions and comments 
helped to improve the paper. 
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